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ABSTRACT 

This is the first of two papers investigating the deprojection and spherical averaging 
of ellipsoidal galaxy clusters. We specifically consider applications to hydrostatic X- 
ray and Sunyaev-Zel'dovich (SZ) studies, though many of the results also apply to 
isotropic dispersion-supported stellar dynamical systems. Here we present analytical 
formulas for galaxy clusters described by a gravitational potential that is a triaxial 
ellipsoid of constant shape and orientation. For this model type we show that the 
mass bias due to spherically averaging X-ray observations is independent of the tem- 
perature profile, and for the special case of a scale-free logarithmic potential, there is 
exactly zero mass bias for any shape, orientation, and temperature profile. The ratio 
of spherically averaged intracluster medium (ICM) pressures obtained from SZ and 
X-ray measurements depends only on the ICM intrinsic shape, projection orientation, 
and ifrji which provides another illustration of how cluster geometry can be recovered 
through a combination of X-ray and SZ measurements. We also demonstrate that Ysz 
and Yx have different biases owing to spherical averaging, which leads to an offset 
in the spherically averaged Ysz — Yx relation. A potentially useful application of the 
analytical formulas presented is to assess the error range of an observable (e.g., mass, 
Ysz) accounting for deviations from assumed spherical symmetry, without having to 
perform the ellipsoidal deprojection explicitly. Finally, for dedicated ellipsoidal stud- 
ies, we also generalize the spherical onion peeling method to the triaxial case for a 
given shape and orientation. 

Key words: X-rays: galaxies: clusters — X-rays: galaxies — dark matter — cosmo- 
logical parameters — cosmology observations 



INTRODUCTION 



Although it is well-known that galaxy clusters are not spherical, spherically averaged measurements of cluster properties are 
standard practice owing to their comparative simplicity and the expectation that non-spherical effect s do not dominate the 
error budget. However, as cosmolo gical measurements with clusters become ever more precise (e.g., IVikhlinin et al.l 1200°] ; 



Allen et al.ll201ll : IPierre et allboilh , the need for precise control of systematic errors also increases. Deviations from spherical 



symmetry, such as intrinsic flattening and substructure, will necessaril y introduce scatter, and possibly biases , into spherically 
averaged global scaling relations used in c osmological studies (e.g., Evrard et al. 19961 : White et al. 20021 : Kravtsov et al. 
2006; IShaw et al.l 120081 iKrause et al.l 120111 ) e.g., the scaling between total mass and average intracluster medium (ICM) 



temperature, quantities that are usually computed interior to a spherical volume specified by a standard fraction of the virial 
radius. 

Unfortunately, owing to their increased complexity and computational expense, non-spherical models are rarely accounted 
for in detail, if at all, in the error budgets of cluster measurements. The modeling of a cluster possessing substructure 
requires sophisticated three-dimensional N-body simulations, while even t o compute the gravitational potential of simple 
ellipsoidal mass distributions involves solving complicated integrals (e.g., Chandrasekhar 1987 : Binnev fc Tremainel 20081 ) 



* E-mail: buote@uci.cdu 



2 D. A. Buote and P. J. Humphrey 



unless app roximate fitting formu las a ppropria t e for nearly spherical systems are employed ( Lee fc Sutol2004 ). The pioneering 
studies by Piffaretti et al, ( 20031 ) and lGavazzi 1 20051 ) dedicated to the problem of assessing errors arising from the assumption 
of spherical symmetry emp loyed ellipsoidal models with some simplifying ass umptions to reduce the total computational 
expense; e.g.. iGavazzi (2005) studied the face-on projections of spheroidal NFW (jNavarro et al.ll 19971 ) mass distributions with 
isothermal ICM. These authors found that quantities obtained from spherical X-ray studies vary typically by < 5% as a result 
of different intrinsic shapes and viewing orientations for a cluster. While not a large effect in absolute terms, errors of a few 
percent may be important for precision cosmological studies and deserve further elucidation. 

The non-spherical shapes of clusters need not be considered merely a nuisance as a source of systematic error, since 
they are interesting to study in their own right. With the aid of large, cosmological N-body simulations signifi cant theoretical 



progress has been made in the understanding of the intrins i c shapes of ACDM cluster dark matter halos (e.g., Jing et al 



2004: 



1995; 



Mohr et al.lll995l; ISplinter et al.lll997l; iThomas et al.|[l99sl; Ijing fc Sutcll2002l; lBullockll2002l; [Suwa et al" book ISpringel et al 



Hopkins et al. 2005; 



Kasun fc Evrard 



2005; 



Lee et al 



2006: IShaw et al. 2006; Bett et al.ll2007l : lGottlober fc Yepe s 2007; Maccio et al.l l2008: Mun oz-Cuartas et al]|201ll ; iRossi et al 



2005; lAllgood et al. 



120061; iLeel I2006T 



Paz et al.1 |200^T|Ho et al 



2011). The axial ratios of dark matter halos in ACDM are found to be sensitive to the value of erg, the normalization of the 
power spectrum of density fluctuations, and are weakly sensitive to the cosmological matter density parameter, fi m . Dark 
matter halos in ACDM become flatter with increasing mass, and cluster-mass halos tend to be prolate/triaxial. Furthermore, 
halos will be more spherical if the dark matter particle i s self-interacting, which provides an additional constraint on particle 
dark matter models, especially on the galaxy scale (e.g., Feng et al. 20091 . 2010l ). 

For a fe w clusters X-ray obse r vations have measured flattened dark matter halos by adopting hydrostatic models 
of the ICM ( Fabricant et al, 1984 ; Buote fc Canizares 19921 . 1996b ). and there is recent evid ence for triaxiality i n the 



dark matter when X-ray d ata are combined with constraints from strong gravitational lensing ( Morandi et al. 201fj| . 2011 



Morandi fc Limousin 2011 



Isothermal hydrostatic models of the ICM applied to the dark matter halos formed in cosmologi- 
cal si mulations pre dict distributions of average X-ray isophotal axial ratios consistent with cluster observations (jWang fc Fan! 
20041 ; iFlores et aklboOTl ; iKawaharalkoid ) . When baryons are included in N-body, hyd rodynamical simulations of cosmic struc- 



ture, it is found that baryon condensation leads to rounder dark matter halos (e.g.. |Pubinskilll994l ; iKazantzidis et al.l 2004; 
Debattista et al. 20081 ). However, the treatment of baryon evolution in cosmological simulations remains a problem for cluster 



cores. For example, baryon cooling in some cosmol ogical simulation s can lead to highly flattened ICM cluster cores that 
disagree with X-ray observations of relaxed clusters (jFang et al. 20091 ) . whereas the global ICM shapes of the model clusters 
agree with X-ray observ ations and also with the shapes of the total gravitat ional potentials of the models as expected for 
hydrostatic equilibrium ( Buote fc Tsai 1995 ; Fang et al]|2009l; [Lau et al.l 2011). IC M shapes have also been measured via the 
thermal Sunyaev-Zel'dovich (SZ) effect l|De Filippis et alj|2005l ; ISavers et al.ll201ll ) which provides another promising avenue 
for probing intrinsic cluster shapes. 

The non-spherical results cited above represent only a small minority of X-ray and SZ cluster studies. To a large extent 
spherical models dominate because they are easier to implement and are substantially more computationally efficient than are 
models which involve a gravitational potential generated by an ellipsoidal mass distribution. Here we investigate a different 
type of model where the potential, rather than the underlying mass distribution, is an ellipsoid of constant shape and 
orientation. These ellipsoidal models lend themselves to straightforward generalizations of simple, analytic spherical models, 
and are therefore just as computationally efficient. Moreover, we show that many cluster quantities derived assuming spherical 
symmetry can be easily interpreted in terms of these ellipsoidal models for a given shape and orientation. 

The paper is organized as follows. In Sj2]we define the ellipsoidal models. We show in fj3]that the relationship between 
the mass profile and potential for these models closely resembles the spherical case, as does that between the mass profile and 
the ICM density and temperature for the case of hydrostatic equilibrium. In [J4]we provide several analytic expressions for 
quantities associated with these ellipsoidal models, in particular relating their deprojected spherical averages to their intrinsic 
ellipsoidal profiles. For practical use of these models with binned data, in [J5] we generalize the traditional spherical onion 
peeling deproje ction method appropriate fo r a series of concentric, triaxial ellipsoidal shells. Our conclusions are presented in 
31 In Paper 2 |Buote fc HumphrevlEnibl ) we perform a detailed investigation of biases and scatter in the measurements of 
global quantities resulting from the spherical averaging of ellipsoidal galaxy clusters. Finally, we mention that the formulas 
we present here (and the results in Paper 2) apply not just for massive clusters but also for groups and massive elliptical 
galaxies with hot gaseous halos. 



2 ELLIPSOIDAL POTENTIALS 

Consider an ellipsoid with principal axes a, b, c and axis ratios, p v = b/a and q v = c/a, satisfying < q v < p v < 1. When the 
ellipsoid is aligned so that a lies along the a;-axis, b along the y-axis, and c along the z-axis, then the ellipsoidal radius a v is 
given by, 
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pi 



(i) 



We define an ellipsoidal gravitational potential (EP) to be an ellipsoid of constant shape (p v ,q v ) and orientation so that $ 
depends only on a v ; i.e., $ = &(a v ). Since $ has a constant shape, EPs best approximate those clusters having nearly constant 
ICM shapes. While it is well-known that individual clusters, both those observed and formed in cosmological simulations, can 
sometimes exhibit large radial shape variations, we note that the ave rage X-ray ellipticity profile for a low-redshift cluster 
sample varies only weakly with radius (see Fig. 8 of Fang et al. 20091 ) . In contrast, the shapes of the isodensity surfaces of 
the underlying mass distribution of an EP generally vary with radius. This can b e a desirable featur e since dark matter halos 
formed in cosmological simulations typically have radially varying shapes (e.g., Jing fe Suto 20021 ; Kazantzidis et al. 2004 ; 



Bailin fc Steinmet jliobll ; Ullgood et aflbood ; IVera-Ciro et alJboilh 



Probably the key advantage of an EP is that a simple, analytic form for $ can be adopted based on a straightforward 
generalization of a spherical potential, enabling much faster computational evaluation than for ellipsoidal mass distributions; 
e.g., see Paper 2 for EP generalizations of the NFW and isothermal j3 models. A disadvantage of an EP is that when the 
potential is sufficiently flattened (the amount depending on how steep is the rad i al potential profile), the mass density can 
become negative in some region (e.g.. iBinnevI Il98ll ; iBlandford fc Kochanekl Il987l ; I Schneider et al.l 1 19921 ; iKassiola fc Kovner 

iej[20! 



19931 ; lEvans! 1994; Bi nnev fc Tremain Requiring a non-negative phase-space distribution function (DF) will further 



restrict the flattening of the matter dis tribution that gen erates the potential if the DF depends only on the energy and one- 
component of the angular momentum ( Evans! 19931 . llQQ^ ) . although it is unclear that the same restrictions would apply to a 
general three-integral DF. Moreover, we show below f i]3.1[) that the mass enclosed within a v has a form analogous to the mass 
enclosed within the spherical radius r. It is therefore as well-behaved as the mass profile of a spherical m odel, indicating that 
EPs should be suitable fo r many applications. Indee d, in our previous X -ray studies of elliptical galaxies ( Buote fc CanizaresI 



1996a ; Buote et al. 20021 ) and a simulated cluster ( Buote fc Tsai 1995j) we found very good agreement between the gravi- 



tating mass profiles inferred using EPs with those obtained using the far more computationally expensive ellipsoidal mass 
distributions. 



3 MASS DISTRIBUTION 
3.1 General Case 

As noted previously, unlike $ itself, the mass density p(x, y, z) of an EP is not an ellipsoid, and its typically complicated form 
must be inferred from direct solution of Poisson's Equation. However, M(< a v ), the mass enclosed within a v , is far simpler 
to compute given $, as we now show. 

Gauss's Law for the mass enclosed within a surface S is, 

Menc = * V$ ■ MS. (2) 

« s 

Now assuming <!> = $(a„) for an EP and taking S to be the ellipsoidal surface defined by a v , the gradient and vector surface 
element take the form, 

V*(a«) = (xl+^j+^k) — and MS = I XQvi = + Vqvj = + k)dxdy. (3) 

\ Pv Qv J a v aa v WaJ-a; 2 -!/ 2 /^ pi ^ a 2 , — x 2 — y 2 /pl ) 

Substituting these expressions into eqn. © and rearranging into two separate surface integrals yields, 



y 2 



M(<a v ) = — / / dxdy + — — - — — / / £=• ~ . dxdy. (4 

4ttG a v da v J J a 2 - x 2 - y 2 jpl 4ttG a v da v J J pi \pi qi J a 2 - x 2 - y 2 /p 2 

Because of the ellipsoidal symmetry we may evaluate the integrals in any quadrant of the x — y plane. Multiplying the result 
obtained from one quadrant by four gives the total surface integral in the positive z-direction, and multiplying by eight gives 
the result for the total ellipsoidal surface. By effecting a change of variable, y = up v \J 'al — x 2 , the first term on the R.H.S. in 
eqn. Q simplifies to, 



2 p v q v d$ 
7iG a v da v 



dx 



2 1 2 

qi I qi 



1 du _a% d$ p v q v ( 2 , 
VT^ 1? ~ G da v 3 I ql ' 



Here the u-integral evaluates to tt/2 so that the a;-integral is elementary. Similarly, making the same change of variable in the 
second term on the R.H.S. in eqn. ((4]) yields, 

2 p^ d$ n _ i_\ n , 2 _ a) r 1 = «| d$ ^ n _ A 

ttG a v da v [pi ql) J v > J jY~^ G da v 3 \p 2 q 2 J w 
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where this time the it-integral evaluates to 7r/4, again leading to an elementary x-integral. Substituting the results from 
equations ([5]) and ([6]) back into eqn. Q gives the result for the enclosed mass, which we state as a theorem: 

Theorem 1. For an EP the mass enclosed within the ellipsoid of radius a v is, 

M(<a v ) = r)(jp v ,q v )^^-, (7) 

where we define the EP shape factor, 
1 1 



, v _ Pvqv 
f]{p v ,q v ) = — — 



1 + 



(8) 



which is of order unity. 



Hence, the relation between M(< a v ) and <&(a v ) is nearly identical in form to the spherical case, M(< r) = (r 2 /G)d$/dr, to 
which it reduces for p v = q v — 1. Here it is worth emphasizing that the axis ratios p v and q v define the shape of <3> and the 
bounding surface for the calculation of M(< a v ), but generally the axis ratios of the corresponding mass density distribution 
are smaller (e.g., see Fig. 1 of Paper 2). 

The spherically averaged mass distribution of an EP does not possess such a simple relationship to the potential. Nev- 
ertheless, it can be computed without resorting to the evaluation of the volume integral of the generally complicated density 
distribution. 

Theorem 2. For an EP the mass enclosed within the sphere of radius r is, 

, s r f d$(a v ) r f M(< a v ) 
M{< r) = — - / a v — f — -dfi = - — r / — i -dll, (9) 

4tG J 4v da v 4nr l (p v ,q v ) J 47T a v 

where dQ = sin 8d9dcf> is the solid angle and the integration proceeds over the entire spherical surface. 
Proof. Again using Gauss's Law (eqn. [2}, but this time for a spherical surface, gives, 

MS = rr 2 dS} = (xi + yj + zk) rdQ, (10) 

for the vector surface element. Taking the dot product with V'l > (a„) (eqn. [3} and integrating over the sphere gives the stated 
result, where the R.H.S. made use of Theorem [T] □ 



3.2 Hydrostatic Equilibrium 



Since a primary goal of our study is to investigate the effect of spherical averaging on the inferred mass distribution from 
X-ray observations, we need to consider the case where the hot intracluster medium (ICM), or "hot gas", is in hydrostatic 
equilibrium. We assume the self-gravity of the gas can be neglec ted ($ ga a $)) w hich is generally a very good approximation 
interior to rsoo where the cluster gas fraction is 
equilibrium, 



10% (e.g.. IPratt et alJboioh . In this case, the equation of hydrostatic 



VPgaa = -p g asV<E>, 



(11) 



where P gas is the thermal pressure and p gas is the density of the ICM, requires that surfaces of constant potential are the 
same as surfaces of constant ICM pressure, density, temperature, and, so long as the metal abundances do not vary over 
these surfaces, it follows that the surfaces of constant X-ra y emissivity also follow the potential ("X-ray Shape Theorem," 
Buote &i Canizaresl 11994 . Il996al ; iBuote &i Humphrevll2011al ). For the special case of an EP model it follows that all ICM 

fsas = Pga.s(a v ), Pgas = Pg as (a v ), and T = T(a v ), so that the hydrostatic equation becomes, 

(12) 



dR 



sas (a v ) _ _ . . d<b(a v ) 

— Pgas \O v J - 



da v aa v 

where we have used the definition of the gradient in eqn. (|3]). Using Theorem [T] and substituting the ideal gas equation of 
state for the pressure, P gas = PgasfcsT/(pw-a), where fcs is Boltzmann's constant, m a is the atomic mass unit, and /x is the 
mean atomic weight of the gas, we obtain the following result. 



Theorem 3. For an EP in hydrostatic equilibrium the mass enclosed within an ellipsoid of radius a v is, 



M(< a v ) = -r)(p v ,q v ) 



a v k B T 



pm a G 



+ 



dlnT 



d In p ga 

din an dlna„ 



(13) 



where rj is given by eqn. (J8j> . 
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This result is a simple generalization of the spherical case where r is replaced by a v and the mass is multiplied by the shape 
factor r/(p v ,q v ). Similarly, solutions of eqn. (1 12[) for p ga s, T, and the entropy are also easily constructed from their spherical 
counte rparts by making the transformation r — > a v and M(< r) — > M(< a v )/r/ (e.g., see equations 10-12 of lBuote fc Humphrey 
2011al ). 



4 DEPROJECTION RELATIONS FOR SPHERICALLY AVERAGED QUANTITIES 

Here we derive analytical expressions for the spherically averaged deprojection of intrinsically ellipsoidal quantities, with a 
particular application to the EPs. We begin by summarizing the general results for the projection of a triaxial ellipsoid. 



4.1 Preliminaries 



We define the orientation of the ellipsoidal system following Binnev 1 1985h . In the reference (x,y,z) coordinate system the 
principal axes of the ellipsoid are aligned with the coordinate directions as described in the definition of a v (eqn. [T]). Define 
a new rotated coordinate system (x',y',z') where the z axis lies along the line-of-sight to the observer from the center of 
the ellipsoid and the x' axis is located in the (a;, y) plane. The two systems are related by starting with their axes aligned, 
rotating th e reference sy stem first by an angle cj> about the z axis and then by an angle 9 about the x' axis. (For the rotation 



matrix, see 



Binnevl[l985l . 



For any quantity that depends only on the ellipsoidal radius a v , Contopoulos! (1956) and Stark (1977) showed that the 
projection of this ellipsoidal quantity yields a (two-dimensional) elliptical distribution of constant shape and orientation that 
depends only on the elliptical coordinate a 3 on the sky; i.e., the projection of the ellipsoidal volume emissivity e(a v ) (i.e., 
luminosity per volume) yields an elliptical surface brightness E(o s ), 



EM = 77 



e(a v )a v da v 



(14) 



where, 



f ■ 2 a I 2 i sax 
j = sm 9 I cos <p + 



pi 



+ 



cos 2 8 

9 > 



(15) 



using the angle definitions described above, and where p v and q v define a v . The elliptical coordinate variable a s — "/ s qaCe is 
proportional to the elliptical radius a (i.e., semi-major axis on the sky), which is of more immediate interest to the observer, 



a — X H 

qi 



(16) 



where q s is the axial ratio of E, and X and Y are sky coordinates aligned with the isophotal major and minor axes respectiveljQ. 
The proportionality factors are given by, 



1 

2? 



(A + C) + A/(A-Cr + B 2 



2 _ (A + C) — V(A-C) 2 + B 2 



(A + C) + v/(A-C) 2 + B 2 
which were derived by StarkI (1977), though also see iBinggelil |l980h and lBinnevI (|1985T ). where 

sin 2 9 



■ 2 , COS 

sm m + 



Pi 



+ 



2 ' 

pi 



l \ l 

cos 9 sin 2(p [ 1 ^ I — , 

Pv J Qv 



pi 



+ COS 



qi ' 



(17) 
(18) 

(19) 
(20) 
(21) 



To convert our notation to that used bv lStarkl jl977f) for the projected quantities, let 7 S — > a, q s — > /3, a — > q, and X 4$ Y. 
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4.2 General Case 

Definition 1. For an elliptical distribution S(a), where a is the elliptical radius with axial ratio q s , we define the "effective 
circular average" (E(7?)} by associating the circular radius with the geometric mean radius: (E(i?)} = E(a), where R = a^/ql- 

The use of such an effective circular average is required for obtaining the analytical relations we describe below. We note 
that for typical models (i.e., those considered in Paper 2) the effective circular average usually very closely approximates a 
formal azimuthal integration of E(a) at fixed R. Only for very flattened models that fall steeply with radius do we find notable 
differences. We emphasize, however, that for our purposes it is only necessary that the observer employs the effective circular 
average in their analysis (as is a very common practice). 

Definition 2. For any quantity intrinsically distributed as an ellipsoid stratified on surfaces of constant ellipsoidal radius 
a v , let (■ • -} d represent the deprojected, spherically averaged version of that quantity defined in the following manner. For 
illustration consider the specific case of e(a v ) and its projection E(a s ) as defined previously. The spherical averaging considered 
here begins by first taking the effective circular average of E(a 3 ) (i.e., Definition [TJ and using the fact that a s oc a) to give the 
radial profile, (E(i?)), where R is the radius on the sky. Then (E(i?)) is deprojected by assuming spherical symmetry, which 
yields the spherically averaged volume emissivity denoted by, (e(r)) d . 

Definition 3. For any quantity that is composed of one or more deprojected spherically averaged quantities as defined in 
Definition^ but itself does not satisfy Definition^ we shall also refer to this composite quantity as a "deprojected spherical 
average" but add a plus sign in the notation to distinguish it from Definition [2] i.e., (• ■ -) d+ . 

Theorem 4. For any ellipsoidal distribution that depends only on the ellipsoidal radius a v , the deprojected spherical average 
(Definition [2} of this distribution is, 

(e(r)) d = ^7 s g s 2 e ( a v), where a v = (-y s ^/q^)r. (22) 

Proof. We start with e(a v ) and use eqn. (|14|l to compute E(a s ). From Definition [T] the effective circular average of E(a s ) is 
obtained by associating the circular radius R with the geometric mean radius of the ellipse of semi-major axis a, R — a^fql. 
Since the ellipsoidal coordinate a 3 is a function of the elliptical radius, a s (a) = 7 s Q s o, it follows that, 

<E(i?)> = E(a,(a)) = E(a s (-p=)) = ^(isVqiR), (23) 

which relates the circular distribution on the L.H.S. to the elliptical distribution on the R.H.S. Since (E(ii)) is circularly 
symmetric, it may be spherically deprojected using the inverse Abel integral relation, 

^ 1 r d ^ R » - dR (24) 

Changing the integration variable from R to a s yields, 

(e(r) = — / — ■ = , (25) 



where a v — {"isy/q7)r. Since eqn. (|14|) is an Abel integral, its inverse is readily obtained (e.g., eqn. 13 of IStarkHl977r ) 



f \ V7 f dE(a s ) da s 
e{a v ) = / — . (26) 

Comparing equations (|25[) and (|26|) gives the desired result. □ 
4.3 X-Ray Emission and Hydrostatic Equilibrium 

Here we consider cluster properties associated with the ICM X-ray emission. We assume that all volume ICM properties depend 
only on ellipsoidal radius a v , which applies exactly for hydrostatic equilibrium in an EP (see i]3.2[) . However, hydrostatic 
equilibrium is strictly required below only for Theorems [7] and and Corollary [2] 

Theorem 5. The deprojected spherical average of the emission-weighted temperature is, 

(T(r)) d+ = T(a v ), where a v = {j s ^)r. (27) 

Proof. Let e be the X-ray emissivity and T the gas temperature. The emission-weighted temperature is denned as the volume 
integral of eT divided by the volume integral of e. Hence, at any radius, 

{(fT\(r)) d (isql f~^\ e{a v )T{a v ) 
(T(r)) d+ = U JZl - V -, i \, = T(a v ), (28) 
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where a v — ('y sy /q^)r, and we made use of Theorem [4] in both the numerator and denominator. Note that the deprojected 
spherical averages in the numerator and denominator each correspond to Definition [2] while the composite result {T(r)) d+ 
corresponds to Definition [3] Similar usages apply below. D 



Theorem 6. The deprojected spherical average of the gas density is, 
(p g as(r)) d+ = (js qi f~ z ^ p gas (a t ,), where a v = ijsy/ql)r. 
Proof. From the definition of the X-ray emissivity, 

d+ _ r < e M> d 



A((T(r)> 



d+1 



(29) 



(30) 



where we have suppressed the metallicity dependence of the plasma emissivity A, which does not affect our arguments provided 
that the metallicity depends only on a v as we are assuming for all ICM properties. Applying Theorem [4] in the numerator 
and Theorem [5] in the denominator of the above equation yields, 



<p gas (r)) d+ = 



yysQsf *) e(a v ) 
A(T(a„)) 



Is q 



V 4 



3 (a v )\(T(a v )) 



A(T(a„)) 



(31) 



where a v — ('Jsy/qs)r, which reduces to the desired result. □ 

Corollary 1. The deprojected spherical averages of the gas pressure and entropy are, 

(P gas (r)) d+ = (tJ^ /"*) PgasM, {S(r)) d+ = (77*97*/*) S(a v ), where a v = (j s <Jqs~)r. (32) 

Proof. These results are immediate consequences of the definitions of each quantity, P gas = p S askBT/(fj,m a ) and S = 
(fcs//im a )Tpg a y 3 , Theorem[S] and TheoremjS] □ 



Theorem 7. The deprojected spherical average of the total mass enclosed within radius r is, 
(M(< r)) d+ = {^y 3 qi r](p v ,q v )j M(< a v ), where a v = {"fs^/ql)r, 
for any temperature profile. 

Proof. Applying hydrostatic equilibrium for a spherically symmetric cluster gives, 



(33) 



<M(< r)) d+ = 



r 2 d 



d + 



r 2 d 



(34) 



<p g as(r)) d + G dr x * Basv 11 p S zs(a v ) G dr 

where a v = ( r y s ^cfJ)r, and the R.H.S. made use of Corollary [1] Theorem[6] and the fact that y s , q s , and / depend only onp„, q v 
and the fixed line-of-sight projection orientation. Changing the variable from r to a v = (7sy / g7)r so that d/dr = ( r y s ^q^)d/da v 
after simplifying, gives, 

— 1 a 2 d 



<M(<r)) d+ = 



1 



~1sJq s 



Pgas (a v ) G da v 



(35) 



By making use of eqn. (|12l) and Theorem[T]the quantity in brackets equals M(< a v )/rj(p v ,q v ), which proves the theorem, for 
any temperature profile T(a v ). □ 

We have chosen to state explicitly that this result for the total mass holds for any temperature profile since we desire to 
emphasize this point below in Theorem [9] Next, however, we consider the gas mass. 



Theorem 8. The deprojected spherical average of the gas mass enclosed within radius r is, 

d+ f ~ ~ 1 \ _1 

(M gas (< r)> = l-Js 2 qs* f 4 Pvq v J M gas (< a v ), where a v = (■Jsy/q^r. 
Proof. The gas mass within a spherical volume of radius r is, 

>gas(a„)47rr 2 dr, 



<M gas (< r)) d+ = / (p gas (r)) d +4^r 2 dr = (7.*^/"*) 



(36) 



(37) 



where a v = {js^/qs)r, and in the R.H.S. we applied Theorem [6] and again (as in Theorem]?} made use of the fact that 7 S , q s 
and / do not depend on r. Changing the integration variable from r to a v = (7 s y / 97)f , and simplifying, gives, 



(M gas (< r)) d + = 



i (a v )4iia 2 da v . 



(38) 
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Since an ellipsoidal volume element is dV = 4np v q v alda v , the integral equals M gas (< a v ) j (p v q v ), which proves the theorem. 
□ 

Corollary 2. The deprojected spherical average of the gas mass fraction enclosed within radius r is, 

d+ / 2 3 i \-l 

(/gas(< r)) = yyiql f i Pvqv/r](pv,qv)\ /gas (< Ov), where a„ = (jax/ql) r. (39) 

Proof. This result follows immediately from the definition of the gas fraction, / gas (< r) — Af gas (< r)/M(< r), Theorem [71 
and Theorem [U □ 

Theorem 9. For the scale-free logarithmic EP, &(a v ) — (GM A /a A ) ln(a„), spherical averaging does not bias the mass profile 
in the sense that, 

(M(< r)) d+ = (M(< r)) truo , (40) 

where (M(< r)) truc is the spherical average of true mass distribution corresponding to $(a„). This result is independent of 
the gas temperature profile. 

Proof. This is a special case of the general problem that is the focus of Paper 2; i.e., we wish to compare spherically averaged 
quantities obtained by an observer to those obtained by theoretical studies. While an observer will measure deprojected 
spherically averaged quantities (Definitions [2] and |3j , the theorist typically spherically averages the true three-dimensional 
distribution directly. From Theorem [2] we obtain the spherical average of the mass distribution generated by the scale-free 
logarithmic EP, <&(a„) = (GMa/ua) ln(ai,), where M A and a A are constants, 

(M(< r)) truc = / av ^±dn = M A — . (41) 
47rG L da v a A 

•/ 47T" 

For comparison, the deprojected spherical average of the mass profile (Theorem [Jj) depends on the mass enclosed within 
ellipsoidal radius a v , which is obtained by inserting the definition of the scale-free $(a„) into eqn. (JjJ of Theorem [TJ 

M(< a v ) = r](p v ,q v )M A — = (laqlvb?v,qv)) M A — , (42) 

OA V / OA 

where the R.H.S. made use of the substitution a v — (js^fqH)r as appropriate for Theorem[jJ Now substituting M(< a v ) into 
Theorem [7J gives, (M(< r)} d+ = M A r/a A = {M(< r)) tTua , independent of the gas temperature profile, which proves the 
theorem. □ 

This theorem complements and extends the result presented in Appendix B of Churazov et al. 1 20081 ). These authors 
consider the bias due to spherical averaging of an isothermal ICM with a scale-free gas density, p gas = /i(S,0)r _Q! , where 
h(8, (f>) is some positive function. They argue that this gas density distribution leads to an inferred total mass that also 
displays no bias due to spherical averaging. Because their model implies a potential, $ oc ln[h(8, <f))r~ a ], equivalent to the 
scale-free potential we employed above in Theorem [9] for a = 1 and h — l/y/J, it is reassuring that the two results each 
predict no bias for the special case of an isothermal ICM. The assumption of the scale-free EP with no a priori restriction on 
the form of the gas density has allowed us to generalize rigorously the zero-bias result for any ICM temperature profile. 



4.4 Sunyaev-Zel'dovich Effect and Related X-ray Quantities 

Now we consider a galaxy cluster also to be observed via the thermal SZ effect, and we continue to assume that all three- 
dimensional ICM properties depend only on the ellipsoidal radius a v . 

Corollary 3. The deprojected spherical average of the ICM electron pressure, P e = uJ^bT, obtained from a measurement 
of the thermal SZ effect is, 

(Pe(r))tz = (jsqlf'^J Pe{a v ), where a v = (ja-s/qTir. (43) 

Proof. The thermal SZ effect is the temperature decrement, ST/T = —2y c in the Rayleigh- Jeans tail of the Cosmic Background 
Radiation (CBR) spectrum due to inverse Compton scattering of CBR photons by energetic ICM electrons. The Compton-y 
parameter is, 

<Tt f „ , i a T 2 f°° P e (a v )a v da v 

Vc= 2 I Psdz = 2^T / / ■> ■> ' ( 44 ) 

m = c2 7ios m - c V7 Ja s Vat - ai 

where or is the Thomson cross section, m e is the electron mass, and the R.H.S. follows from the condition that P e = P e (a v ) 
(i.e., eqn. U4[) . By associating y c (a s ) with E(a s ) and or P e (a„) / (m c c 2 ) with e(a„) the result follows immediately from Theorem 

® n 
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The geometrical factor j s qi /~ 5 in eqn. (|43[1 is the square of the corresponding factor for the ICM pressure obtained 
from X-ray studies (Corollary [l}, indicating that spherical averaging has a stronger impact on the ICM pressure inferred from 
SZ studies. Since this difference in geometrical factors can be exploited to measure the intrinsic shape and orientation of a 
cluster, we state it formally. 



Corollary 4. 

studies is, 



The ratio of the deprojected spherical averages of the ICM electron pressures inferred from SZ and X-ray 



<P S z,x(r)) d+ = 



(Pe(r)) d S z 



(45) 



Proof. This result follows immediately from Corollaries [T] and O where we have used the ICM electron pressure from X-rays, 
<ftM)x + = ((2 + M)/5)(P gas (r)) d +. □ 

We remark that (Psz,x(r)} d+ is, in fact, constant with radius and does not overtly display any dependence on the distance 
to the cluster, and hence the Hubble Constant, Ho- However, Px oc p gas oc y/e, where e is the volume emissivity. To obtain 
physical units for the emissivity requires converting the observed X-ray flux to a luminosity density, the net result of which 
is that Px is inversely proportional to the square root of the cluster distance, leading to (Psz,x(>")) d+ K 1/V-Hb- 

The possibility of uncovering the intrinsic shape of the cluster ICM by combining X - ray and SZ measurements has been 
recognized for over ten years (e.g., lz"aroubi et al.lll99Sl; IFqx fc Penlbooj iReblinskvlboOOl : iLee fc Suto1l2004l ; be Filippis et all 
20051 ; IPuchwein fc Bartelmannlkood ; ISerenof 120071 ; lMahda v^^^ChMig| [2011]) Thi s promising t echnique has already provided 



interesting constraints on cluster shapes for many clusters (|De Filippis et al 



20051 ; ISereno et alj|2006h using isothermal triaxial 



/3 mode ls for the ICM, and more recently models with a radially varying temperature profile have been applied to the cluster 
A 1689 ( Sereno et al. 201 lh . Corollary [4] defines a particular approach to this problem that has some attractive characteristics. 
First, the ratio of spherically averaged pressures is valid for any temperature profile T(a v ). Second, the relationship does not 
assume a particular ICM radial density profile (e.g., /3 model) and, in principle, can be deprojected using the traditional 
spherical onion peeling procedure. Hence, studies can be conducted entirely in the context of spherical symmetry to obtain 
(Psz,x(r)) d+ and then, supplemented with a measurement of the average ICM axial ratio on the sky (q s ), can be used to 
constrain the geometrical factor jg g* / 4 . For the general triaxial ellipsoid this factor depends on the intrinsic shape via the 
axial ratios p v and q v and the orientation (6,(j>). For spheroids only a single axial ratio and inclination angle are required. 
We now consider the quantity, 



Ysz 



(It 



PedV, 



(46) 



where V is the volume and Da(z) is the angular diameter distance to the cluster. Since this quantity equals the integral of 
y c over solid angl e in the limit of a small angle subtended on the sky, it is usually referred to as the "integrated Compton-y 
parameter" (e.g., White e t al. 2002). We will take V to be a large spherical or ellipsoidal region centered on the cluster. 



Theorem 10. The deprojected spherical average of the integrated Compton-y parameter is, 
{Y S z(< r)} d+ = (jsqsf zpvqv^l Y sz (< a v ), where a v = (^ a y/q^)r. 



Proof. Starting from the definition of Ysz, we have for a spherical volume, 



(Ysz« r)) 



d+ 



(JT 



D\(z) m e c 2 
(7.?.*/~ 4 ) 



(P e (r))i z 47vr 2 dr 



(Tt 



P e (a^)47rr dr, 



(47) 



(48) 



(49) 



D\{z) m e c 2 Jo 

where we have used Corollary [3] and the fact that j s , q s , and / depend only on p v , q v , and the fixed line-of-sight projection 
orientation. Effecting a change of variable within the integral from r to a„ = (j s y/qH)r yields, 

1 1 



{Ysz(< r)) c 



= Hqsf- 



(Tt 

D 2 A (z) m e c 2 Jlt 



Y sz (< a v ) 



P e (a v )4na v da v 



p v q v 



(50) 
(51) 



where the last step used the definition of an ellipsoidal volume element, dV = 4np v q v a^da v , which proves the theorem. □ 

Cosmological simulations predict a strong correlation between Ysz and cluster mass (e.g., IWhite et al]l2002r i. which is a 
direct result of the gas pressure probing the depth of the cluster potential wel l. However, becau se simulations do not perfectly 
match observations of cluster ICM (e.g., isophotal flattening of cool cores, iFang et alJboogh . and since it is preferable to 
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use an independent method to constrain cosmology, there is interest in using X-ray observations of cluster mass to calibrate 
Ysz independently. For high-quality X-ray data of clusters where it is possible to measure accurately the spatially resolved 
gas density and temper ature profiles, direct calculation of the ICM pressure profile is to be preferred for comparison to 
Ys ? llArnaud et aL 201C). For lower quality data it is necessary to rely on scaling relations, such as Yx, a quantity advocated 
bv lKravtsov et al as a mass proxy for cosmological studies, 

Y x {<r) = M gas {<r)T x (<r), (52) 

where M gas (< r) is the gas mass and Tx(< r) is the emission- weighted temperature enclosed within the spherical volume of 
radius r, and typically a radius rsoo is adopted. Since Yx is related to the integrated gas pressure profile, it should be closely 
related to Ysz. Indeed, recent SZ studies find a strong correlation between Ysz(< rsoo) and Yk(< rsoo) and between Ysz(< rsoo) 
and M gas (< rso o), each corre l ation having similar intrinsic scatter (e.g., Andersson et al. 2O10l ; Planck Collaboration et al 



201 lh - see also lFabian et al.l |201lh . 



Before addressing the spherical average of Yx, we consider Tx. 

Theorem 11. The deprojected spherical average of the emission- weighted temperature integrated over the spherical volume 
of radius r is, 

(T x (< r)) d+ = T x (< a v ), where a v = ^ s ^q7)r. (53) 

Proof. Beginning with the definition of (Tx(< r)) d+ as the deprojected integrated emission-weighted temperature, we have, 

\Z{((-T){r)) d A-Kr 2 dr f r e(a v )T(a v )A-Kr 2 dr 

(T x (< r)) d+ = J "V A " = Jo V ^ V ; , (54) 

Jo {^(r)} d 4nr 2 dr J Q e(a v )Anr 2 dr 



where we made use of Theorem 0] in both the numerator and denominator and the fact that the factor (7s y/Qs/f) depends 
only on p v , q v , and the fixed line-of-sight projection orientation. Effecting a change of variable within the integrals from r to 
a-v = ('Ja^ql)r yields, 

(T x (< r)> = a „ . r: — j-, ■ (55) 

J Q €(a v )4-Kaida v 

Multiplying the numerator and denominator of the R.H.S. by p v q v yields the emission-weighted temperature integrated over 
the volume of ellipsoidal radius a v = {'isyfqTjr, which is the desired result. □ 

We remark that Theorem [11] reduces to Theorem [5] for the special case of a small radial volume element associated with 
a finite radius r. 

Corollary 5. The deprojected spherical average of Yx(< r) is, 

(Y x (< r)) d+ = (jlq^f^PvqJ) Y x (< a v ), where a v = ("fsjq^r. (56) 
Proof. This result is an immediate consequence of the definition of Yx(< r) (eqn. [52)) and Theorems [8l and 1111 □ 
Corollary 6. The ratio of the deprojected spherical averages of Ysz and Yx is, 

(Y sz ,x(< r)) d+ = ^y^^C = «?/"*) y x Z (<£) = T l ) Y S z,x(< a.), where a v = ( 7 ^)r. (57) 

Proof. This result follows immediately from Theorem [10] and Corollary [5] □ 



4.5 Connection to Stellar Dynamics 

Many of the results we have presented can be applied either directly, or with minor modification, to a relaxed, dispersion- 
supported collisionless stellar system with an isotropic velocity dispersion tensor. Such a system obeys the equation of hy- 
drostatic equilibrium where the stellar density p s ta.rs(a, v ) replaces the ICM density and the square of the velocity dispersion 
a(a v ) 2 replaces the gas temperature. (Here it is assumed the stars, like the gas, are merely a tracer of the gravitational 
potential, which is a good approximation for galaxy clusters, and also for elliptical galaxies well outside of the stellar half-light 
radius.) Consequently, (p s tars(r)) d obeys Theorem[4l and (cr(< r) 2 ) d+ obeys Theorem[5] Similarly, the deprojected spherically 
averaged mass inferred from the stellar dynamics also obeys Theorem [7] where in the proof one replaces the gas density with 
Pstars(rau) an d the pressure with Pstars(au)o"(a„) 2 . Finally, the deprojected spherically averaged stellar mass profile behaves as 
(Ysz(< r)) d+ (Theorem I10p . because the deprojected stellar mass density behaves as (P e (r))g Z (Corollary [3} • 
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5 PROJECTION OF ELLIPSOIDAL SHELLS AND ONION PEELING DEPROJECTION 

To treat the case of binned observational data, such as the one-dimensional surface brightness profile of a cluster, here we 
describe the projection and deprojection of a system of concentric, similar triaxial ellipsoidal shells relevant for the study of 
EPs. For an ellipsoidal shell defined between aj, n and a° ut with constant emissivity, e(a^, o° ut ), throughout the shell, equation 
(|14|l becomes, 

EKVrVs) = 2£(a '"_ a ° Ut) [ V(aS ut ) 2 - a! - V(< n ) 2 -ail , (58) 
V / L J 

where a s < a™*, and the second term in brackets is set to zero if a s > ajf. This equation projects a three-dimensional 
ellipsoidal shell onto a two-dimensional elliptical surface brightness that depends only on the elliptical coordinate a s . We 
desire the luminosity integrated over an elliptical annulus defined between semi-major axes, ot m and a° ut : 

/■os out 

L(a™, a° ut ; a m , a out ) = / £(a"\ h™'; a a )2-Kq a ada (59) 

J Q,in 

/ in out\Trint/ in out in out\ f nr\\ 

= ^{a v ,a v )V \a v ,a v ;a ,a ), (60) 

where, 

r rint / in out in out \ 47T 

V (a v ,a v ;a , a ) = -^-p v qv x 

0/ out\2 / in\2l3/ 2 ["/ out\2 / out\2l3/ 2 . \f in\2 / out\21 3 /2 ["/ in\2 / in\2~| 3 / 2 \ 

(a v ) - {-Ysqsa ) J - [(a v ) - (7 3 g a a ) J + [(o„ ) - (~f 3 qsa ) J - [(o„ ) - ("/sq s a ) J I . (61) 



If any terms in equation (|61|l have negative arguments, they must be set to zero. Note when viewed "edge-on" (9 = <j> = 90°), 
i.e., down the intermediate principal axis, ellipsoids always have q s — q v and "/ a q s = 1. For the special case of an oblate 
spheroid (p v — 1) viewed at arbitrary inclination, q s — q v yfj and 7 s q s = 1. 

Let an ellipsoid be partitioned into a series of concentric, similar ellipsoidal shells, a v ^ < a Vl i < a„,2 < ■ • ■ < a.v,N- Define a 
corresponding set of concentric, similar, elliptical annuli such that, ao < u\ < 012 < ■ ■ ■ < Qn, where ao = i»,o, ot\ = a v ,2, ■ ■ ■■ 
For this case, we may represent the projection of the three-dimensional ellipsoidal shell (aj-i,aj) onto the two-dimensional 
elliptical annulus (aj_i,Qj) by, 

V]f = V lnt (a v j-i,a v j;ai-i,ai). (62) 

That is, each shell and annulus is labeled by the index of its outer boundary. The contribution of shell j to the luminosity of 
annulus i is, Li = ejVj^, where tj = e(a v ,j-i,civ,j) is the constant emissivity within the shell. We obtain the total luminosity 
projected into annulus i by summing the contributions from all shells j > i, 



N 



tiVjT, (63) 



so that the surface brightness is, 

T 1 N 

^ = (2 2 S = 7 2 2 S E ' W'"- (64) 

Hence, we have shown that the projection matrix for spherical shells (e.g., equation B12 of lGastaldello et~al 120071 ) is generalized 



to the case of ellipsoidal symmetry via the following mapping: three-dimensional radius, r — > a v ; two-dimensional radius, 
R — > a s = 7 s g s a; and 47r/3 — > p v q v 4ir/3. Moreover, by separating the first term from the others in the summation of 
eqn. (f63|) . 

N 

Li = eiVjV 1 * + ^ ej'^'i 1 *, (65) 
and then solving for the emissivity in shell i, 



we arrive at a generalization of the "onion peeling" deprojection method ( Fabian et al. 1981 : Kriss et al.ll 19831 ) appropriate for 



triaxial ellipsoids (p v , q v ) with any orientation (6, cj>). That is, the emissivity in shell i is obtained by taking the total luminosity 
observed in annulus i on the sky, subtracting from it the luminosity contributions projected from shells at larger radii (j > i), 
and then dividing by V/" representing the volume of intersection between the ellipsoidal shell i and the ellipsoidal cylindrical 
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shell defined by the base of elliptical annulus i and an infinite height. Practical implementation of eqn. (|66[1 requires assuming 
values for 6, <f>, p v , and q v and measuring the value of q s on the sky. 

While the above has focused on the example of the emissivity projecting into the surface brightness, generalization to 
other quantities such as the emission-weighted temperature and projected temperature map ((T)j) is straightforward; i.e., 

N 

(T^i-^e/Z^f. (67) 

3=i 



6 CONCLUSIONS 

This is the first of two papers investigating the deprojection and spherical averaging of ellipsoidal galaxy clusters (and 
massive elliptical galaxies). We specifically consider applications to X-ray and SZ studies, though many of the results also 
apply to isotropic dispersion-supported stellar dynamical systems. A major disadvantage of working with ellipsoidal systems, 
as opposed to spherical systems, is that they generally involve numerical evaluation of computationally expensive integrals. 
Here we present analytical formulas for galaxy clusters described by a gravitational potential that is a triaxial ellipsoid of 
constant shape and orientation; i.e., an "ellipsoidal potential" (EP), which depends only on ellipsoidal radius, $ = $(a„). 
While the mass density is itself not ellipsoidal for these models, and it can take unphysical values in the vicinity of the 
minor axis when the flattening is too large, we demonstrate that the total mass enclosed within the ellipsoidal radius a v is 
proportional to d<&/da v , and is therefore well-behaved for any smooth &(a v ), making it useful for many purposes. 

We show that for hydrostatic X-ray studies of EPs the relationship between the enclosed total mass, ICM temperature, 
and ICM density has the same form (up to a proportionality factor) as the spherical case where the spherical radius r is 
replaced by a v . Using this result, along with the general result we derive for the spherical deprojection of any ellipsoid of 
constant shape and orientation, we show that the mass bias due to spherically averaging X-ray observations is independent 
of the temperature profile. For the special case of a scale-free logarithmic EP ($ oc lna„) there is exactly zero bias for any 
shape, orientation, and temperature profile. The ratio of spherically averaged ICM pressures obtained from SZ and X-ray 
measurements depends only on the intrinsic shape and projection orientation of the EP, as well as Ho, which provides another 
illustration of how cluster geometry can be recovered through a combination of X-ray and SZ measurements, with the key 
advantage that the pressures are measured in the context of spherical symmetry without (in principle) having to specify a 
parametric form for the radial profile. We also demonstrate that Ysz and Yx have different biases as a result of spherical 
averaging, which lead to an offset in the spherically averaged Ysz — Yx relation. Paper 2 explores in more detail the biases 
and scatter arising from spherical averaging, in particular using the observationally and cosmologically motivated NFW mass 
profile, and also considers the more widely investigated, and computationally expensive, class of potentials where the mass 
density, rather than the potential itself, is an ellipsoid of constant shape and orientation. 

A potentially useful application of these analytical formulas is to assess the error range on an observable accounting for 
deviations from assumed spherical symmetry without having to perform the ellipsoidal deprojection explicitly. That is, an 
X-ray observer can, as is standard, analyze a cluster assuming spherical s ymmetry and obtain deprojected ICM temperature 
and density profiles using the spherical onion peeling procedure (see fj5j Fabian et al. 1981 ; Kriss et al. 19831 ) to construct, 



e.g., the observed spherically averaged mass profile, (M(< r)) d+ . With the aid of Theorem [7] this can be converted into the 
true mass profile M(< a v ) for an assumed three-dimensional shape and viewing orientation. Then by using Theorem [2] the 
true spherically averaged mass profile {M(< r)) truo can be constructed. By adopting priors for the shape and orientation and 
marginalizing over them, the full range of {M(< r)) true owing to intrinsic ellipsoidal geometry can be computed. 

Finally, for dedicated ellipsoidal studies, we also generalize the spherical onion peeling method to the triaxial case for 
a given shape and orientation. The formulas presented for ellipsoidal shells may also be of use for ellipsoidal projections in 
numerical work. 
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